The Ji — J2 Heisenberg model at and close to its 2; = 4 quantum critical point 
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We study the frustrated Ji — J2 Heisenberg model with ferromagnetic nearest neighbor coupling 
Ji < and antiferromagnetic next-nearest neighbor coupling J2 > at and close to the z; = 4 
quantum critical point (QCP) at J1/J2 = —4. The Ji — J2 model plays an important role for 
recently synthesized chain cuprates as well as in supersymmetric Yang Mills theories. We study 
the thermodynamic properties using field theory, a modified spin-wave theory as well as numerical 
density-matrix renormalization group calculations. Furthermore, we compare with results for the 
classical model obtained by analytical methods and Monte-Carlo simulations. As one of our main 
results we present numerical evidence that the susceptibility at the QCP seems to diverge with 
temperature T as x ~ 2 jjj ^jjg quantum case in contrast to the classical model where x ~ 

PACS numbers: 03.67.-a, 11.25.Hf, 71.10.Pm, 75.10.Jm 



I. INTRODUCTION 



The frustrated one-dnnensional (ID) s — 1/2 Heisen- 
berg model 



H - Ji ^ SjSj+i 



(1) 



with nearest neighbor coupling Ji and next-nearest 
neighbor coupling J2 is the minimal model to describe 
magnetism in a number of cuprate chain compounds. It 
can be viewed as a ladder with coupling J2 along the 
legs and a zigzag rung coupling Ji as shown in Fig. 1 (a) . 
Recently, there has been renewed interest in this model 
with ferromagnetic coupling Ji < and antiferromag- 
netic coupling J2 > propelled by the discovery of multi- 
ferroic behavior in edge sharing spin chains. The phase 
diagram of this model as a function of a = Ji/ J2 has been 
studied by a combination of field theoretical and numer- 
ical methods. ^"^'^ It is shown schematically in Fig. 1(b). 
At a = the system consists of two decoupled, critical 
antiferromagnetic Heisenberg chains. By bosonization it 
has been found that a small coupling, |Ji| <C 1, leads 
to an exponentially small gap A cx exp(— const |ci;|).®'^'^^ 
On the antiferromagnetic side, a > 0, the gapped phase 
exists up to a critical point (QCP2 in Fig. 1(b)) at 
a « 4.15"'^^ where the system enters a critical gapless 
phase. At the so-called Majumdar- Ghosh (MG) point, 
(5 = 2, the ground state is known exactly and consists 
of decoupled dimers. Dimerization is indeed present for 
the whole gapped phase < a < 4.15 while short-ranged 
incommensurate (SRI) spin correlations have only been 
found for < a < 2.^'^'^ On the ferromagnetic side, 
(5 < 0, a phase with incommensurate spin-spin correla- 
tions is followed by a ferromagnetic phase. The transition 
occurs at a = — 4 (QCPl in Fig. 1(b)), both in the quan- 
tum as well as in the classical model. ^^'^^ Whereas the 
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FIG. 1: (a) The Ji — J2 chain viewed as a ladder with zigzag 
coupling, (b) Phase diagram of the Ji — J2 chain as a function 
of J1/J2 with J2 > 0. 



incommensurate ("spiral") correlations are long- ranged 
in the classical model, these correlations are predicted 
to be short-ranged in the quantum model. However, 
the gap is expected to be exponentially small and no nu- 
merical evidence for this gap has been found yet.^'""^*'^^ 
At the critical point QCPl, the ferromagnetic state and 
states of resonance valence bond (RVB) character are 
degenerate. ^° In fact, all degenerate ground states at this 
point can be explicitly constructed. ^"'^'^^ It turns out that 
there exists a unique ground state for a fixed total spin 
5tot and fixed z-component of the total spin, S^^^. The 
ground state with S^^^ = is, in particular, a uniformly 
distributed RVB state obtained as a superposition of all 
possible states where sites are grouped in singlet pairs. 

In this paper we will study the thermodynamic proper- 
ties of the Ji — J2 model near the quantum critical point 
QGPl. There are two reasons why this model is of cur- 
rent interest: On the one hand, the recently studied com- 
pound Li2ZrCu04 has been shown to be well described 
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by the Ji — J2 model with a frustration parameter a 
putting the system into the spiral phase but rather close 
to the critical point QCPl.^'^^ By chemical or external 
presssure it might be possible to tune this or a related 
system across the phase transition. On the other hand, it 
has been shown using the anti de Sitter/conformal field 
theory (ADS/CFT) correspondence that a deep connec- 
tion between spin chains and string theory exists. ^^'^^ In 
A/" = 4 super Yang Mills theories the dilatation opera- 
tor in two- loop order can be represented as the S = 1/2 
spin chain, Eq. (1), with parameters fixed by the Yang 
Mills coupling constant. ^'"''^^ In the relevant parameter 
regime, both couplings Ji and J2 are ferromagnetic in 
this case. Interestingly, however, the second order con- 
tribution taken separately has J1/J2 = —4 although with 
J2 < 0, i.e., there is no frustration. We will see in Sec. II 
that for this specific ratio — irrespective of the sign of J2 — 
certain terms in the effective field theory will cancel ex- 
actly. 

Our paper is organized as follows. In Sec. II we present 
a field theoretical description of the model in the ferro- 
magnetic phase. Based on this field theoretical model 
we will also discuss the properties of the critical point 
QCPl when approached from the ferromagnetic side. In 
Sec. Ill we investigate the thermodynamics of the clas- 
sical model. We test the analytical results for the low- 
temperature properties obtained in Sec. II by comparing 
with Monte Carlo (MC) simulations. In Sec. IV A analyt- 
ical results for the quantum model based on a modified 
spin-wave theory (MSWT) are obtained. We compare 
the MSWT predictions with the field theory and with 
numerical data obtained by the density-matrix renor- 
malization group algorithm applied to transfer matrices 
(TMRG) in Sec. IV B. A summary and conclusions are 
presented in Sec. V. 



II. FIELD THEORY AND SCALING 
ARGUMENTS 

We consider the case a < 0. Using spin coherent 
statcs^^ the Hamiltonian (1) can be mapped onto a non- 
linear sigma model with Euclidean action 

Se = -isJ^Swz[n{r)] + SH (2) 

r 

with (we set ft = Ajb = 1) 

Sh = dTy2{Jis'^n{r,T)n{r + ao,T) 
Jo J. 

+ J2s'^n{r,T)n{r + 2ao,T)}. (3) 

Here n^(r, r) = 1 is a unit vector, s the spin quantum 
number and ao the lattice constant. Swz['n{r)] is a topo- 
logical (Berry) term giving a phase which is determined 
geometrically by the cap bounded by the trajectory n{T). 
Without the topological term we have a classical action. 



Parametrizing the unit vector in terms of angle variables 
and demanding that the action is stationary we can eas- 
ily find the classical ground state. This leads to the well- 
known result that the ground state is ferromagnetic for 
a < — 4 and a spiral with pitch angle = arccos(|a|/4) 
for a > —4. 

Up to a constant we can replace n(r, T)n{r + ao,T) — >■ 
— (n(r, r) — n(r -I- ao,T))^/2 and similarly for the next- 
nearest neighbor term. In the continuum limit wc can 
then expand the action in terms of the lattice constant 
ao and obtain in leading orders 



S„ = - 



J2S^ao 



+ 
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J2s'^al 
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dr 



dr{'i + a){drnf 



I'dT [ dr{16 + a){dSf , (4) 
Jo Jo 

with L = Nao where N is the number of lattice sites. 

It is instructive to briefly discuss the planar case 
where n is restricted to the x-y plane. In this case 
we can parametrize the unit vector by a single angle, 
n = (cos (j), sin (p, 0), leading to 



S 



H 



J2S^ao 
2 

24 







/3 fL 

d,T / dr{4: + a){dr(pf 
Jo 

dT I dr(16 + a) [{dl(l)f + {drcj))' 




(5) 



For a < —4 we see that the action is minimized by 
dr4> = d^.(fi = 0, i.e., the ground state is ferromagnetic. 
For a > — 4 the system can gain energy by forming a 
spin spiral. Right at the transition point, the first line of 
Eq. (5) vanishes. The dispersion of the elementary exci- 
tations (first term in the second line of Eq. (5)) therefore 
becomes quartic at the critical point QCPl, ojk ^ k*- 
The critical theory therefore has a dynamical critical ex- 
ponent 2 = 4 whereas in the ferromagnetic phase the 
dispersion is quadratic {z = 2). The dynamical critical 
exponent relates the scaling of energy to and length L, 
CO L~^. For the free energy / = — ^ In Z it follows that 
/ ~ — T'3/2 ^Yie ferromagnetic phase and / ^ — T^/* at 
the critical point. The same scaling relations apply for 
the inner energy. For the specific heat c = —Td^f/dT^ it 
follows c ^ T^/^ in the ferromagnetic phase and c ^ T^/^ 
at the critical point. These scaling relations will stay 
valid also in the general case described by Eq. (4) and 
depend only on the dimension of the dynamical critical 
exponent. 

Next, we consider the magnetic susceptibility in the 
ferromagnetic phase. The operator in the second line of 
Eq. (4) is then irrelevant and can be ignored. The par- 
tition function, including a magnetic field h, to leading 
order is then given by 



-I 



Dn exp ■ 



dr 



{drnf - hMon^ 
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,(6) 
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with the spin stiffness ps = — s^ao(Ji + 4J2) and Mq = 
s/ao- fn the second line wo have rcscaled r' = Tr/ps 
and introduced a new parameter g = hM^pg/T"^ . We are 
always interested in TL/ p^ ^ 1, i.e., in systems at tem- 
peratures T much larger than the finite size gap ~ 1/L. 
If g is the only parameter of the theory, then we expect 
a universal scaling for the magnetization, M = Mq^{()), 
where ^ g + 0{g^) is a universal scaling function. 
For the susceptibility it follows that x ^ M^ps/T^. 

Following Ref. 28 one can even go one step further 
and calculate the scaling function ^{g) explicitly. To do 
so it is important to realize that Eq. (6) is nothing but 
the imaginary time path integral of a quantum particle 
moving on a sphere. The corresponding Hamiltonian is 
then given by H = 1? j1 — gn^ where L is the angular 
momentum operator. The scaling function can now be 
obtained by calculating the eigenspectrum of this Hamil- 
tonian leading to ^{g) = + 0{g^).'^^ If the scaling 
hypothesis is valid, we expect the susceptibility at low 
temperatures on the ferromagnetic side of the transition 
( Ji + 4 J2 < 0) to be given by 



X = 



3 r2 



2.s4 Ji + 4 J2 



3 



rp2 



' 3T2 



a 



(7) 

This relation has been found first in Ref. 17 based on 
an analysis of numerical data. The low-temperature be- 
havior is therefore the same as for the nearest-neighbor 
ferromagnetic Heisenberg model^^ but with a rescaled 
spin stiff'ness ps- At the critical point QCPl we have 
Ps ^ signalling the formation of spiral correlations. In 
our treatment we have ignored the Berry phase term. In 
analogy to the simple ferromagnetic model, we expect 
that this term in the ferromagnetic phase does not play 
any role for the low-temperature physics and thus the 
low-temperature thermodynamic properties of the quan- 
tum and the classical J1-J2 model are the same. 

Let us now consider the field theory, Eq. (4), at the 
critical point , still ignoring the topological term in Eq. 2. 
The partition function is then given by 



Z — i Dn exp 



1 



dr 



Ps {dlnf-hM^n'- 



dr' 



{dl,nf-gn^ (8) 



with Ps = |Ji|s^og/8 and g = hMopV^ /T'^/^. The sus- 
ceptibility therefore scales as 



n^2;;l/3 8/3 



l^ll 



2-4/3 



2 



(9) 



where C is a constant. Two of us have shown that one 
can again go one step further by considering Z as the 
path integral of a quantum anharmonic oscillator. The 
eigenvalues of this Hamiltonian can then be calculated 
numerically. In full analogy to the ferromagnetic case 



discussed before, the proportionality factor can therefore 
be determined and is given numerically by C « 2.14. 

The result (9) is expected to be the exact low temper- 
ature susceptibility for the classical model at a = —4. 
However, for the quantum model the theory is above the 
upper critical dimension, d + z = 5, with d being the 
dimension of the system. Spin-wave interaction terms 
might therefore yield ultraviolet (UV) divergencies so 
that the result for the susceptibility might not only de- 
pend on the parameter g but also on a UV cutoff. "^^ In 
this case the scaling hypothesis would be violated and 
the formula (9) not applicable for the quantum model. 
Furthermore, the topological term which we have ne- 
glected throughout, is likely to play an important role at 
QCPl. Here the ferromagnetic state is degenerate with 
RVB states^" which do not exist for the classical model. 
Within the nonlinear sigma model description one might 
expect that part of this difference is encoded in a nontriv- 
ial topological term. In the following, we will first check 
the analytical predictions for the classical model before 
analyzing a possible violation of Eq. (9) for the quantum 
model. 



III. THE CLASSICAL MODEL 

The classical nearest-neighbor ferromagnetic Heisen- 
berg model has been solved by Fisher [31]. For the clas- 
sical model (1) with antiferromagnetic Ji and J2, Harada 
and Mikeska [32] have shown that thermodynamic quan- 
tities can be expressed in terms of eigenvalues of transfer 
matrices which follow from integral equations. The scal- 
ing X ~ T^^/"^ at the Lifshitz point, a — —4, has already 
been discussed in Ref. [33] . Recently, also the thermody- 
namics for general Ji < and J2 > has been studied 
in more detail. 2^'^^"^'' 

In the classical case, the results derived in the previous 
section by field theory methods — both for the ferromag- 
netic phase and the critical point — should be valid. Here 
we concentrate on providing numerical evidence that the 
parameter-free formulas for the magnetic susceptibility, 
Eq. (7) and (9), respectively, are correct. The numerical 
data are obtained using Monte Carlo simulations from 
the ALPS package^^ with a cluster update and a system 
size N = 10000. 

For the ferromagnetic phase, Eq. (7) predicts that at 
low temperatures all data for x/[s'^{Ji+4:J2)] should col- 
lapse onto a single universal curve. In Fig. 2 Monte Carlo 
results for various a are compared with the analytical 
formula. The data collapse onto the analytical curve is 
perfect over temperatures of several orders of magnitude. 
We note that the closer a is to the critical point QCPl 
the lower the temperatures are where the universal scal- 
ing sets in. 

Similarly, we can also check the formula for the critical 
point, a = —4, see Fig. 3. The numerical data do confirm 
the analytical result, however, we note that temperatures 
T/l Ji| < 0.01 are required. 
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FIG. 2: (Color online) Universal data collapse, x/i^'^i^J'^ ~ 
4J2)] = 2/(3r^), at low temperatures for the classical model. 
The line denotes the analytical result (7), and the symbols 
MC results for a — —5 (circles), 3 = — 10 (squares), and 
a = —20 (diamonds). 
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FIG. 3: (Color online) Analytical formula (9) for the low- 
temperature susceptibility at a = —4 (line) compared to MC 
data (circles). 



IV. THE QUANTUM MODEL 



A. Modified spin wave theory 

To calculate the thermodynamic properties of the spin- 
s ferromagnetic Heisenberg chain, Takahashi [39,40] in- 
troduced a modified spin-wave theory. The spin oper- 
ators are represented by bosonic operators as in reg- 
ular spin-wave theory. In addition, the constraint of 
vanishing magnetization at finite temperatures posed by 
the Mermin- Wagner theorem is implemented in a sim- 
ple way by adding an effective magnetic field which acts 
as a Lagrange multiplier for the magnetization. MSWT 
has also been used successfully to describe bound- 
ary contributions in the open ferromagnetic Heisenberg 
chain*^ as well as the thermodynamics in the dimer- 
ized ferromagnet.'*^"'** Moreover it has been shown that 
the classical ferromagnetic chain is well-described by 
MSWT.^O''*^'^^ Here we will apply the same method to 
the Hamiltonian (1) with general spin s. 

We expect that MSWT can be apphed for a < — 4 
where the ground state is ferromagnetic. We will also 
use this approximation for the QCP at a = —4, how- 
ever, here the validity of MSWT is questionable because 
of the degeneracy of the ferromagnetic ground state with 
RVB states. A further discussion will be presented in 
Sec. IV B, based on a comparison with numerical data, 
and in the conclusions, Sec. V. In the MSWT approxi- 
mation the Hamiltonian (I) is represented as (in the fol- 
lowing we set the lattice constant oq = 1) 



H = N{Ji + J2)s^ 



k 



(10) 



where a[J' is a bosonic annihilation (creation) opera- 
tor with [afe,a^,] ~ 5k,k'- The dispersion relation is 
given by Uk = 2s(| Ji| [1 - cosfc] + J2 [cos(2fc) - 1]). The 
additional constraint of vanishing magnetization reads 
s = '^f^ rik where n^. = (exp[ajj./T -I- u] — 1)^"'" is 
the Bose function including the effective magnetic field, 
V = h/T. 

First, we study the case a < —4. For small temper- 
atures only spin-wave excitations with small momenta 
contribute and the dispersion can be approximated as 
<^fe — |>^i|s(l ~ 4/|a|)fc^. The constraint can now be 
solved explicitly by expanding in the reduced temper- 
ature t. This leads to 



The quantum s — 1/2 nearest-neighbor (J2 — 0) fer- 
romagnetic Heisenberg chain is exactly solvable. Ther- 
modynamic properties, in particular the susceptibility, 
have been calculated using the thermodynamic Bethe 
ansatz.'^^''"^ The results have been shown to be in excel- 
lent agreement with those obtained by a modified spin- 
wave theory (MSWT).'*° In the following we extend the 
MSWT approach to the ferromagnetic phase of the J1-J2 
model and to the quantum critical point. We then test 
the obtained analytical results by comparing with numer- 
ical data obtained by the transfer matrix renormalization 



group. 
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t = 




2s 



\JMi-4/\a\y 



(11) 



Here Ci^) is the Riemann zeta-function. We note that 
the expression for v is the same as for the simple 
ferromagnet,'*^ only the definition of the reduced tem- 
perature is modified. The expansion (11) is valid if 
\/i/2s <C 1, i.e., the temperature range where this re- 
sult is applicable shrinks the closer we get to the critical 
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point. The free energy in spin- wave approximation is 
given by 



/ = (Ji + J2)s' - T 



^^ln(l + n,) 



(12) 



At small temperatures we can expand the second term 
in V and obtain 



/ = ( + J2)s^ - T 



C(3/2) 

2 V TT As 



(13) 



In agreement with the scaling relations derived in Sec. II 
we find as leading temperature dependence / ~ —T^l^. 
From C = —Td'^f/dT'^ the leading temperature depen- 
dence of the specific heat can be obtained. 
The susceptibility in MSWT is given by 



1 



^nfc(nfe -f 1) 



(14) 



leading to the low-temperature expansion 



X 



3 r 
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(15) 



The leading temperature dependence found in MSWT 
therefore agrees exactly with formula (7) found by gen- 
eral scaling arguments. 

Next, we consider the critical point a ~ —4. For small 
momenta the dispersion now reads ujk — \ J\\sk'^lA and 
thus the dispersion changes from quadratic to quartic 
which will have consequences for the temperature scaling 
of thermodynamic quantities. For the constraint we find 
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f34s4 
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V2C(l/4)p/4 
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with a new reduced temperature t = 4r/(|Ji|s). Using 
relation (12) we obtain for the free energy 



/ = (Ji + J2)s2 _ T 



r(l/4)C(5/4) -i/, 
47r 



t 

Ms 



1/3- 



(17) 



plus higher order terms. The scaling of the leading term, 
/ ~ — y5/4^ jg again consistent with the scaling argu- 
ments in Sec. II. Finally, we can calculate the suscepti- 
bility using Eq. (14) and obtain 
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(18) 

We note that the scaling of the leading term is the same 
as in Eq. (9). However, the numerical prefactor is not 
the same. For the case s = 1/2 we find, in particu- 
lar, from (18) that x = 2-^°^^T-^^^ « 0.0992 T-^/s. 52 
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FIG. 4: (Color online) Numerical data (squares) for the 
singlet-triplet gap A^t for chains of even lengths = 
12, 14, • • • , 36. The solid line is a fit A,t = ao + ai (iV"^)"" 
with ao = -4.8 ■ 10"^ ai = 594.1, and 02 = 0.9994. 



In contrast, Eq. (9) — which does give the correct low- 
temperature behavior of the classical model, see Fig. 3 — 
yields x ~ 0.1685 T"^/"^. First of all, this does suggest 
that at the QCP the classical and the quantum model 
no longer show the same thermodynamic properties at 
low temperatures. However, similarly to the scaling ap- 
proach used in Sec. II one might also question the foun- 
dations of MSWT for the QCP aUtogether. At the QCP 
the ground state is no longer a simple ferromagnet sug- 
gesting that also the excitations are no longer described 
by simple spin waves. Furthermore, UV divergencies in 
contributions from spin wave interaction terms might be 
expected because we are now dealing with a theory above 
the upper critical dimension. An independent test of the 
scaling and the MSWT approach at the QCP can only 
be obtained by unbiased numerical calculations. Such 
calculations will be presented in the next subsection. 



B. Numerical results 

The quantum critical point QCPl at a = —4 is charac- 
terized by a level crossing of a singlet and a fully polarized 
state. Right at the critical point the singlet-triplet gap 
Ast is therefore expected to vanish. This is confirmed by 
the Lanczos calculations for finite size chains shown in 
Fig. 4. We find, in particular, that the singlet-triplet gap 
vanishes as Agt ^ N~^. 

The thermodynamics of the Ji — J2 model has been 
studied previously by a Green's function method^^ and 
by TMRG.'^* In the second approach the one-dimensional 
quantum model is mapped onto a two-dimensional clas- 
sical model with the help of a Trotter-Suzuki decompo- 
sition. It is then possible to express the partition func- 
tion in terms of a transfer matrix for the classical model 
with the free energy depending only on the largest eigen- 
value of this transfer matrix. The transfer matrix is ex- 



6 



H 



11 

M 

U 100 




. * e 




'L — 




If 10 

\ \ 

\W 




\ V 




\ "f 


1 


\ > 0.01 0.1 

\\- T/IJ,I 


I.I. 





0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
T/IJjl 



-0.25 




-0.35 



FIG. 5: (Color online) Susceptibility for a = -20. The 
symbols denote the TMRG data, the solid line a solution of 
(14) where the Lagrange parameter v is determined by solv- 
ing the non-linear equations numerically, and the dashed line 
the leading low-temperature asymptotics (15) extracted an- 
alytically. The inset shows the low-temperature region on a 
logarithmic scale. 



tended in imaginary time direction — corresponding to a 
successive lowering of the temperature — with the help of 
a density-matrix renormalization group algorithm. For 
details concerning the algorithm the reader is referred to 
Refs. [41-44,50,51]. 

Here we want to use the TMRG algorithm to test in 
how far the analytical predictions from the previous sec- 
tion hold for the s = 1/2 case. In Fig. 5 numerical 
data for the susceptibility are compared to MSWT for 
a = —20. If we numerically solve the non-linear equation 
for the Lagrange parameter v then the MSWT prediction 
is in excellent agreement with the numerical data up to 
temperatures of order t/s ^ 1. The formula (14) makes 
use of a representation of the susceptibility in terms of the 
spin-spin correlation function (SiSj) which also includes 
terms quartic in the bosonic operators. In this case, the 
constraint fortunately makes it possible to obtain a fi- 
nal expression which is still only bilinear in the bosonic 
operators. The result presented in Fig. 5 therefore goes 
beyond linear spin-wave theory. The formula for the free 
energy, Eq. (12), is, however, a linear spin- wave expres- 
sion. Here it is not possible to include the quartic terms 
without further approximations because the constraint 
alone is not sufficient to obtain a final expression which 
is only bilinear in the bosonic operators. The MSWT 
results for the free energy are therefore only valid at very 
low temperatures as can be seen by the comparison with 
numerical data, Fig. 6. A similar comparison for the crit- 
ical point is shown in Fig. 7. There is good quantitative 
agreement at temperatures t/s < 1, i.e. T/\Ji\ < 0.0625 
between the numerics and the fully self-consistent solu- 
tion of the MSWT equations. 

For the susceptibility the situation is expected to be 
more complex. The scaling hypothesis used to derive 



FIG. 6: (Color online) The free energy for a = —20, —10, —5 
(from bottom to top). The symbols denote the TMRG 
data, the solid lines the MSWT result (12) using a fully self- 
consistent solution for the Lagrange parameter v. 
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FIG. 7: (Color online) Free energy for the critical point a = 
—4. The symbols denote the TMRG data, the dashed line 
the leading terms in the low-temperature asymptotics (17) 
obtained by MSWT. The solid line represents the fully self- 
consistent solution of Eq. (12). 



Eq. (9) is questionable because at the critical point we 
are above the upper critical dimension. Indeed, we have 
already seen that the predictions from MSWT deviate 
from formula (9) which we have confirmed to be the cor- 
rect result for the classical model. This is contrary to 
the ferromagnetic regime where the MSWT results coin- 
cide at low temperatures with the solution of the clas- 
sical model. A comparison with numerical data. Fig. 8, 
indicates nevertheless an apparently good quantitative 
agreement with MSWT up to temperatures T ^ \Ji\. A 
closer inspection of the low-temperature asymptotics (see 
inset (b) of Fig. 8), however, shows that MSWT is not 
fully consistent with the numerics. Furthermore, inset 
(c) shows that classical and quantum model no longer 
share the same low-temperature properties. If we fit the 
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FIG. 8: (Color online) Susceptibility for the critical point 
a = —4. (a) and (b): The symbols denote the TMRG data, 
the solid line a solution of (14) where the Lagrange parameter 
V is determined self-consistently, and the dashed line the low- 
temperature asymptotics (18). The inset (b) shows the low- 
temperature region on a logarithmic scale, (c) Comparison 
between numerical data for the quantum model (small circles) 
and the classical model (squares). 



0.05 0.245 1.195 

0.025 0.259 1.185 

0.01 0.247 1.193 

0.0075 0.244 1.195 

0.005 0.240 1.198 

0.004 0.237 1.201 



TABLE L Parameters obtained by fitting x{T) in Fig. 8 to 
X{T) = A ■ T"' in the interval T/\Ji\ £ [0.003, Tmax]. 

numerically obtained x{T) for the quantum model to a 
simple power law and vary the fit region we obtain the 
values summarized in table I. This seems to indicate 
that the exponent might actually be smaller than 4/3 
and therefore different from the exponent in the classical 
model. However, using the TMRG algorithm we are not 
able to reach temperatures for the quantum model which 
are as low as those obtainable for the classical model us- 
ing MC simulations. We therefore can not completely 
rule out that the temperatures are just not low enough 
to observe the T~^/^ power law predicted by MSWT. 



V. CONCLUSIONS 

In summary, we have studied the frustrated J1-J2 
Heisenberg chain at and near to its z = 4 critical point 
at a = J1/J2 = —4. By developing a field theory we 
have discussed how the system is driven from a ferro- 
magnetic state to a state with incommensurate (spiral) 
spin-spin correlations. Based on this analysis, the clas- 



sical and the quantum model are expected to show the 
same low-energy properties in the ferromagnetic phase. 

From scaling arguments we obtained, in particular, 
that the susceptibility diverges x ~ {Ji + 4J2)/r^ in 
the ferromagnetic phase with a known proportionality 
constant. We could reproduce this result using an alter- 
native analytical approach based on a modified spin wave 
theory. Furthermore, we have verified the analytical pre- 
diction using Monte Carlo simulations for the classical 
model and the transfer matrix rcnormalization group for 
the quantum model and have found excellent agreement. 

Right at the critical point, J1/J2 = —4, simple dimen- 
sional analysis allowed us to predict the scaling of the free 
energy and specific heat with temperature. In particular, 
we found that the free energy scales as / - T^/". Our 
modified spin wave theory calculations have confirmed 
this scaling and the obtained parameter-free results have 
been shown to be in good agreement with numerical data 
for the quantum model at low temperatures. Scaling 
arguments can also be used to obtain a parameter-free 
formula for the low-temperature behavior of the suscep- 
tibility, which — according to this formula — diverges as 
X Ymia the field theory analysis even the pref- 

actor can be obtained and we have shown, using Monte- 
Carlo simulations, that this formula is indeed correct for 
the classical model. 

The most interesting problem is the temperature 
dependence of the susceptibility at the critical point 
Ji/ J2 — —4 for the quantum model. The modified spin 
wave approach also yields a r~'*/3 divergence, however, 
the prefactor is different from the one obtained from field 
theory. The numerical data, furthermore, seem to in- 
dicate that even the exponent might deviate from 4/3. 
From fits of our numerical data at the lowest accessible 
temperatures we have obtained x ~ T~^'^ . 

There are a number of possible reasons for this devi- 
ation: The simplest explanation is that we do not have 
numerical data for low enough temperatures to observe 
the true scaling behavior. However, there are good rea- 
sons to believe that the observed deviation has physical 
reasons. An inspection of the quartic term describing the 
interaction of spin waves within spin wave theory shows 
that this term is ultraviolet divergent. Such divergences 
are expected because d -t- 2: = 5 is larger than the upper 
critical dimension. In such a case the ultraviolet proper- 
ties of the theory can affect the critical behavior. 

Another problem is the treatment of the Berry phase 
term which we have ignored in our field theory analysis. 
On the basis of a bosonization approach,^' where the 
system is considered starting from the decoupling point 
a — Q, the spiral phase of the quantum model has been 
found to be gapped. Since we do not expect an additional 
phase transition, a gap should exist all the way to ci — >■ 
—4. In the field theory approach this gap in the quantum 
model must be related to the Berry phase term. This 
term might therefore also be important for the physical 
properties right at the transition point. This expectation 
seems to be consistent both with the known degeneracy 
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of the ferromagnetic and resonance valence bond states 
at this point and our numerical results which show that 
the low-temperature properties of the classical and the 
quantum model are different. 

The results attained here might be relevant and should 
be compared to future experiments on edge-sharing 
cuprate chain compounds. 
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